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1.  Background 


When  testing  the  performance  of  vehicles  against  under-body  threats  such  as  those  used  in  Iraq 
and  Afghanistan,  the  response  of  a  vehicle’s  structure  is  measured,  in  part,  with  accelerometers. 
Specifically,  measurements  are  taken  from  the  seats  and  floor  as  these  locations  are  in  direct 
contact  with  the  occupant  and  thus  give  some  insight  into  the  loading  condition  experienced. 

This  structural-response  data  can  be  combined  with  measurements  taken  from  anthropomorphic 
test  devices  placed  at  crew  positions  to  provide  valuable  information  regarding  the  interaction 
between  vehicle’s  components  and  occupants.  These  data  can  be  used  to  characterize  the 
environment  that  the  mounted  occupant  is  exposed  to  and  inform  changes  to  the  design  of 
vehicular  components  to  enhance  crew  survivability. 

Drawing  insights  from  accelerometer  data  is  not  necessarily  intuitive  as  the  quantitative  and 
qualitative  characteristics  of  the  data  are  largely  influenced  by  their  measured  frequency  content. 
This  content,  not  readily  discernible  through  visual  inspection,  is  influenced  by  sample  rate  and 
the  need  to  prevent  aliasing  by  implementing  hardware  filters. 

This  point  is  illustrated  with  a  simplified  example  of  a  test  article  exposed  to  a  blast  event.  The 
response  of  the  fictitious  test  article,  shown  in  Fig.  1,  is  composed  of  superimposed  accelerations 
with  magnitudes  of  1,000  m/s  and  10,000  m/s  and  frequencies  of  50  Hz  and  10  kHz, 
respectively. 


Fig.  1  The  contrived  response  of  a  test  article  exposed  to  a  blast  event 
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Environments  with  accelerations  between  300  g  to  1,000,000  g  and  including  mixed  frequency 
content  between  10  Hz  and  1  MHz  are  classified  as  “ballistic  shock”  by  MIL-STD-810G,  a 
military  standard  used  during  Live  Fire  Test  and  Evaluation  (LFT&E)  tests  for  underbody  blast/ 
Though  fictional,  this  example  meets  these  criteria  and  can  therefore  be  considered  a  ballistic- 
shock  event.  The  measurement  techniques  used  to  characterize  it  should  therefore  include 
instrumentation  with  a  “flat  frequency  response  to  at  least  100  kHz”,  a  digitizing  rate  “at  least  2.5 
times”  the  anti-alias  (AA)  filter  frequency,  and  AA  filtering  that  provides  “an  attenuation  of  50 
dB  or  greater  for  frequencies  that  will  fold  back  into  the  passband”.^  The  “measured” 
accelerations  are  compared  in  Fig.  2,  where  they  are  sampled  at  100  kHz  and  20  kHz  with  AA 
filters  of  40  kHz  and  8  kHz,  respectively,  each  simulated  as  digital  4-pole  Butterworth  filters. 


X  1  o'*  1 00/40  kHz  Sample  Rate/AA  Filter 

1.5  I - ^ - , - , - ^ ^ - 


0  0.001  0.002  0.003  0.004  0.005  0.006  0.007  0.008  0.009  0.01 

Time  (seconds) 


20/8  kHz  Sample  Rate//\A  Filter 


Fig.  2  Acceleration  measured  with  100  kHz  and  20  kHz  sample  rates 

Not  surprisingly,  the  8  kHz  AA  filter  required  by  MIL-STD-810G  with  the  20  kHz  sample  rate 
removed  the  contribution  of  the  10  kHz  acceleration  while  the  40  kHz  AA  filter  required  with  the 
100  kHz  sample  rate  did  not.  Perhaps  more  importantly  from  an  assessment  standpoint,  the 
restriction  of  data  bandwidth  to  8  kHz  reduced  the  peak  acceleration  by  more  than  90%  as 
compared  to  applying  the  40  kHz  filter.  Since  it  is  possible  for  the  bandwidth  of  LFT&E  data  to 
vary,  this  example  underscores  the  potentially  misleading  nature  of  acceleration  magnitude  as  a 
metric  when  analyzing  structural  data  and  the  value  of  seeking  alternatives. 

For  many  researchers  and  analysts,  velocity  is  an  intuitive  engineering  parameter.  As  the  integral 
of  acceleration,  it  is  also  less  susceptible  to  the  presence  or  absence  of  relatively  high-frequency 
content  and,  therefore,  is  an  attractive  alternative  to  acceleration-based  metrics.  This  is  easily 
illustrated  by  turning  back  to  the  simple  example.  Velocities  for  the  “actual”  and  “measured” 


2 


responses  are  given  in  Fig.  3.  The  expanded  view  from  9.5  to  10  milliseconds  shows  that  the 
reduction  in  maximum  velocity  associated  with  the  lowest  sample  rate  is  approximately  5%. 


Impact  of  Sampling  on  Velocity 


Fig.  3  Velocity  from  the  actual  acceleration  data  (black),  100  kHz 
sampled  data  (red),  and  20  kHz  sampled  data  (blue) 

Not  only  is  the  assessed  structural  velocity  less  sensitive  to  the  nature  of  the  data  and  how  they 
were  captured,  the  use  of  peak  velocity  as  a  metric  for  vehicle  response  is  useful  from  an 

2  3 

occupant-injury  perspective  due  to  the  known  velocity-dependent  nature  of  human  tissue.  ’ 
However,  since  tissue  has  the  characteristic  of  dissipating  energy  over  time,  it  is  also  useful  to 
know  how  quickly  the  assessed  peak  velocity  is  achieved.  The  ratio  of  the  peak  velocity  to  the 
time-to-peak  velocity  provides  a  constant  pseudo-acceleration,  albeit  one  much  less  susceptible 
to  the  effects  of  sampling  and  filtering. 

To  benefit  from  any  metric,  it  is  imperative  that  it  be  assessed  in  a  uniform  manner  whether 
analyzing  large  amounts  of  data  or  making  comparisons  between  test  events.  So  it  is,  too,  for  the 
calculation  of  peak  velocity  and  time-to-peak  velocity.  To  achieve  the  goal  of  calculating  these 
parameters  uniformly,  a  codified  and  automated  process  should  exist  that  is  applicable  to  as 
many  valid  structural  responses  as  possible.  One  such  process  is  described  in  the  next  section 
with  an  example  and  provided  in  the  form  of  a  MATLAB  script  in  the  Appendix. 
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2.  Methodology 


2.1  Calculation  Method 

The  following  steps  detail  a  simple  method  for  calculating  peak  velocity  and  time-to-peak 
velocity  from  an  acceleration  record.  The  first  step  is  to  apply  baseline  shifts  to  the  acceleration 
data.  Following  integration,  the  appropriate  peak  velocity  is  identified.  The  slope  of  the  velocity 
rise  is  then  calculated  and  used  to  determine  time-to-peak  velocity.  Issues  such  as  rotation  of  the 
test  article,  time  misalignment,  aliasing,  clipping,  and  other  data-acquisition  issues  are  not 
addressed  here  but  should  be  avoided  and,  if  possible,  reconciled  prior  to  the  implementation  of 
this  method. 

2.1.1  Step  1  -  Baseline  Shifts 

Two  distinct  baseline  adjustments  are  made  to  the  acceleration  record  prior  to  integration.  In  the 
first,  the  mean  pre-trigger  acceleration  is  subtracted  from  the  entire  data  record  to  remove  any 
systemic  baseline  shift.  If  the  test  article  comes  to  rest  prior  to  the  end  of  the  data  record,  a 
second  baseline  adjustment  is  applied  in  which  the  average  post-trigger  acceleration  is  subtracted 
from  the  post-trigger  data  record  in  order  to  remove  any  residual  velocity."^  If  a  non-zero  velocity 
is  present  in  the  test  article  at  the  end  of  the  data  record,  the  second  baseline  adjustment  should 
not  be  applied.  As  Fig.  4  shows,  these  shifts  do  not  noticeably  alter  the  acceleration  record  but 
provide  a  meaningful  correction  to  the  velocity  trace.  This  correction  is  small  at  early  times  in 
the  velocity  record  but  more  substantial  at  later  times  due  to  the  accumulation  of  error  when 
integrating  the  acceleration  data. 
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step  1  -  Baseline  Shifts 


Fig.  4  Acceleration  and  velocity  from  the  original  (red)  and  baseline 
shifted  (black)  data 


2.1.2  Step  2  -  Locate  the  Peak 

Following  integration  of  the  acceleration  record,  the  peak  velocity  within  a  time  frame  of  interest 
is  identified.  However,  this  point  does  not  always  represent  the  end  to  the  initial  rise  in  velocity; 
for  example,  it  could  be  a  secondary  peak  that  occurs  much  later  in  time  than  the  initial  velocity 
rise,  as  shown  in  Fig.  5.  Choosing  the  secondary  peak  as  the  peak  velocity  would  enlarge  the 
time-to-peak  velocity  and  therefore  suggest  a  milder  loading  on  the  occupant  than  what  actually 
occurred.  This  is  resolved  by  first  finding  all  local  maxima  in  the  velocity  record  that  have  a 
magnitude  greater  than  90%  of  the  absolute  peak.  The  local  maximum  just  prior  to  the  absolute 
peak  is  then  selected  as  the  new  peak  velocity.  This  search  is  repeated  with  the  new  peak  velocity 
from  the  previous  iteration  until  no  local  maxima  exist  that  meet  the  90%  threshold.  It  should  be 
noted  that  90%  was  generally  found  to  be  the  largest  threshold  that  yielded  the  same  local 
maximum  as  analysts  chose  for  a  majority  of  the  LFT&E  cases  examined  during  the  construction 
of  this  method. 
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Step  2  -  Locate  the  Peak 


-1  I - ^ ^ ^ ^ ^ ^ - 1 

0  5  10  15  20  25  30  35 

Time  (milliseconds) 


Fig.  5  The  absolute  (red)  and  final  (blue)  peak  velocities 

2.1.3  Step  3  -  Calculate  Velocity  Slope 

Once  the  final  peak  velocity  has  been  identified,  an  assumption  of  eonstant  aeceleration  is  used 
to  define  a  velocity  slope.  Three  loeations  on  the  velocity  trace  between  the  start  of  the  time 
frame  of  interest  and  the  peak  velocity  must  be  found.  They  are  the  first  data  points  that  have 
velocity  magnitudes  exeeeding  5%,  20%,  and  95%  of  the  final  peak  velocity.  The  first  and  last 
points  are  the  start  and  end  points,  respeetively,  of  the  velocity  slope.  If  the  velocity  record 
exhibits  a  monotonie  rise  between  the  5%  and  20%  points  (no  negative  derivatives  are  present), 
it  is  appropriate  to  set  the  start  of  the  velocity  slope  relative  to  the  difference  between  the  final 
peak  velocity  and  zero  as  shown  in  Fig.  6.  If  there  is  at  least  one  negative  derivative  in  the  5%  to 
20%  time  window  of  the  velocity  record,  the  starting  point  of  the  slope  must  be  recalculated. 
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Step  3  -  Calculate  Velocity  Slope  (No  Negative  Derivatives) 


0  5  10  15  20  25  30  35 

Time  (milliseconds) 


Fig.  6  Velocity  slope  from  5%  (blue)  to  95%  (red)  of  the  final  peak 
velocity 

To  recalculate  the  starting  point  of  the  velocity  slope,  the  data  point  immediately  following  the 
last  negative  derivative  in  the  window  is  identified  and  used  to  replace  zero  as  the  baseline  for 
the  5%  calculation.  In  other  words,  the  velocity  slope  now  travels  from  5%  of  the  difference 
between  the  final  peak  velocity  and  the  point  after  the  last  negative  derivative  to  95%  of  the 
difference  between  the  final  peak  velocity  and  zero  (as  shown  in  Fig.  7). 
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step  3  -  Calculate  Velocity  Slope  (Negative  Derivatives) 


Fig.  7  Velocity  slope  from  5%  (blue)  of  the  difference  between  the 
final  peak  velocity  and  the  point  after  the  last  negative 
derivative  (green)  to  95%  (red)  of  the  difference  between  the 
final  peak  velocity  and  zero 

2.1.4  Step  4  -  Calculate  Time-to-Peak 

Finding  the  start  and  end  times  of  the  velocity  approximation  is  accomplished  by  extending  the 
slope  line  (Fig.  8).  The  start  time  is  the  intercept  between  the  slope  line  and  zero  velocity,  and 
the  end  time  is  the  intercept  between  the  slope  line  and  peak  velocity.  Time-to-peak  is  simply  the 
difference  between  the  end  and  start  times. 
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step  4  -  Calculate  Time-to-Peak 


Fig.  8  The  fall  approximation  of  the  rise  in  velocity  formed  by 

extending  (blue)  the  velocity  slope  (red)  developed  in  Step  3 


3.  Conclusion 


The  simple  method  explained  in  this  report  was  developed  to  be  applicable  to  a  large  number  of 
structural- velocity  traces  produced  during  LFT&E  with  a  broad  range  of  characteristics.  It 
approximates  the  velocity  of  a  structural  component  from  the  constituent  acceleration  record  in 
order  to  determine  peak  velocity  and  the  time  taken  to  achieve  peak  velocity.  Together,  these 
metrics  are  useful  when  assessing  the  response  of  a  vehicle  exposed  to  an  underbody-blast  event 
and,  consequently,  characterizing  the  environment  of  mounted  vehicle  occupants. 
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Appendix.  MATLAB  Script  for  Calculating  Peak  and  Time-to-Peak  Velocities 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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%Velocity  and  Time-to-peak  Script 
%Robert  Spink,  26  June  2014 


%Input  variables  are  as  follows: 

%time  -  vector  of  time  values 
%vel  -  vector  of  velocity  values 

%search_start  -  start  time  for  the  data  region  of  interest 
%ssearch_stop  -  stop  time  for  the  data  region  of  interest 

%FIND  START  INDEX,  STOP  INDEX,  START  THRESHOLD,  MIDDLE  THRESHOLD,  AND  STOP  THRESHOLD 

start_index=f ind (time>search_start, 1 ) ; 

stop_index=f ind (time>search_stop, 1 ) -1 ; 

start_threshold= . 05; 

mid_threshold= . 2 ; 

stop_threshold= . 95; 

peak_compare= . 90; 


%FIND  ALL  LOCAL  PEAKS 


derivative=dif f (vel (start_index : stop_index) ) ; 
peak_vector=zeros (stop_index-start_index, 1) ; 
for  k=2 : length (peak_vector) ; 

if  (  (derivative  (k-1)  <0&derivative  (k)  >=0)  &vel  (k)  <0)  |  (  (derivative  (k-1) ... 
>0&derivative (k) <=0) &vel (k) >0)  ; 

peak_vector (k) =vel (k+start_index-l) ; 

end 


end 


%PROCESS  DATA 

[ vel_100 , vel_100_index] =max (vel (start_index : stop_index) ) ; 
vel_100_index=vel_100_index+start_index-l; 

recurs ive_peak=length (find (peak_vector ( 1 : vel_100_index-start_index+l-l ) >peak_compare*vel_100 ) ) ; 
while  recursive_peak>0 ; 

vel_l 00_index=f  ind  (peak_vector  ( 1 :  vel_100_index-start_index+l-l ) ... 

>peak_compare*vel_100 , 1 , ' last ' ) +start_index-l ; 
vel_100=vel (vel_100_index) ; 

recurs ive_peak=length  (find  (peak_vector  ( 1 :  vel_100_index-start_index+l-l ) ... 
>peak_compare*vel_100) ) ; 

end 

time_100=time (vel_100_index) ; 

vel_95_index=f ind ( vel ( start_index : vel_100_index-l ) <vel_100*stop_threshold, 1 , ' last ' ) +start_index 
vel_95=vel (vel_95_index) ; 
time_95=time (vel_95_index) ; 

vel_20_index=f ind (vel (start_index : vel_95_index-l ) <vel_100*mid_threshold, 1 , ' last ' ) +start_index; 
vel_20=vel (vel_20_index) ; 
time_20=time (vel_20_index) ; 

vel_5_index=f ind ( vel ( start_index : vel_20_index-l ) <vel_100*start_threshold, 1 , ' last ' ) +start_index; 
vel_5=vel (vel_5_index) ; 
time_5=time (vel_5_index) ; 

if  size (find (derivative ( vel_5_index-start_index+l : vel_20_index-start_index+l ) <0 ) ) >0 ; 

last_inf  lection_index=f  ind  (derivative  ( vel_5_index-start_index+l :  vel_20_index... 
-start_index+l ) <0,1, 'last') +vel_5_index; 
last_inf lection_vel=vel (last_inflection_index) ; 

vel_5_index=f  ind  ( vel  ( start_index :  vel_100_index)  <  ( last_inf  lection_vel+  ( start_threshold*... 

( vel_100-last_inf lection_vel ))),!, 'last') +start_index; 
vel_5=vel (vel_5_index) ; 
time_5=time (vel_5_index) ; 

end 

slope=  (vel_95-vel_5)  /  (time_95-time_5)  ; 
peak_time= ( (vel_100-vel_95) /slope) +time_95; 
start_time=time_95- ( (vel_95-0) /slope) ; 
t2p=peak_time-start_time; 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


AA 

anti-alias 

dB 

decibel 

g 

acceleration  due  to  earth’s  gravity 

Hz 

hertz 

kHz 

kilohertz 

LFT&E 

Live  Fire  Test  and  Evaluation 

m/s^ 

meter  per  second  squared 
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RDRL  SLB  E 
C  BARKER 
B  BENESCH 
D  HO  WEE 
M  MAHAEEEY 
H  MILLER 
J  O’BRUBA 
K  WEAVER 
RDRL  SLB  W 
W  MERMAGEN 
S  SNEAD 

2  DIR  USTARDEC 
(PDE)  RDTA  RTI  GSS 
C  EOSTER 
H  PIETSCH 

1  OSD 
(PDE)  DOTE 

N  BROCKHOEE 

1  USAEC 
(PDE)  TEAE  SV 

S  JANKIEWICZ 

1  USMC 
(PDE)  MCSC 

A  PURTELL 

1  USMC 
(PDE)  MCCDC 
R  CROCE 

1  DOD 
(PDE)  IDA 

J  WALBERT 

18  DIR  US  ARE 

(PDE)  RDRL  DPW 
R  COATES 
P  EROUNEELKER 
T  HOLDREN 
R  SPINK 
M  TEGTMEYER 
JTICE 
RDRL  SEE  L 
L  MERKLE 
RDRL  SEE 
N  SHUKLA 
RDRL  SLB 
B  BOWEN 


1  DEEENSE  TECHNICAL 
(PDE)  INEORMATION  CTR 

DTIC  OCA 

2  DIRECTOR 

(PDE)  US  ARMY  RESEARCH  LAB 
RDRL  CIO  LL 

IMAL  HRA  MAIL  &  RECORDS  MGMT 

1  GOVT  PRINTG  OEC 
(PDE)  AMALHOTRA 
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